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We investigate in the equal-time formalism the derivation and truncation 
^ ■ of infinite hierarchies of equations of motion for the energy moments of the 

covariant Wigner function. From these hierarchies we then extract kinetic 
equations for the physical distribution functions which are related to low- 
order energy moments, and show how to determine the higher order moments 
in terms of these lowest order ones. We apply the general formalism to scalar 
and spinor QED with classical background fields and compare with the results 

derived from the three-dimensional Wigner transformation method. 
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I. INTRODUCTION 



Transport theory [|TJ based on the Wigner operator is extensively used to describe the 
formation and evolution of highly excited nuclear matter produced in relativistic heavy ion 
collisions. The Wigner operator can be defined in 4-dimensional or 3-dimensional []7]|| 
momentum space, which we denote by VV(x,p) and W(x, p), respectively Correspondingly, 
there are two formulations for the phase-space structure of any field. Either of these two 
formulations has its advantages and disadvantages. Besides its manifest Lorentz covariance 
(which is very useful from a technical point of view), another characteristic feature of the 
4-dimensional formulation for QCD @-|5| and QED || is that the quadratic kinetic equation 
can be split up naturally into a transport and a constraint equation. The complementarity of 
these two ingredients is essential for a physical understanding of quantum kinetic theory |§ . 
In the classical limit, these two equations reduce to the Vlasov and mass-shell equations, 
respectively. The main advantage of the 3-dimensional approach f7||| is that it is easier 
to set up as an initial value problem: one can directly compute the initial value of the 
Wigner operator from the corresponding field operators at the same time. In the covariant 
frame this is not possible since the covariant Wigner operator is defined as a 4-dimensional 
Wigner transform of the density matrix and thus includes an integration over time. Hence 
in this approach the initial condition for the Wigner operator at very early times must be 
constructed phenomenologically. Some true quantum problems like pair production fH| in 
a strong external field have thus so far been solved only in the 3-dimensional (or equal-time) 
formulation 0JT1]. 

One way |7|] to obtain equal-time kinetic equations which parallels the procedure in the 
covariant formulation is to Wigner transform the equation of motion for the equal-time 
density operator g(x,y). For spinor QED this procedure results in the BGR equations |7j 
for the equal-time Wigner functions. In Ref. Jl2|| we suggested a different derivation which 
is based on taking the energy average of the covariant kinetic equations in the 4-dimensional 
formulation. It exploits the fact that the equal-time Wigner function is the energy average 
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(i.e. zeroth order energy moment) of the covariant one. With this method we showed 
for spinor QED that the direct energy average of the covariant kinetic equations leads, in 
addition to the BGR transport equations for the spinor components of the equal-time Wigner 
function, also to a second group of constraint equations which couple the equal-time Wigner 
function to the first order energy moment of the covariant one. In the classical (H — > 0) limit, 
these additional equations provide essential constraints on the equal-time Wigner function 
and allow to reduce the number of independent distribution functions by a factor of two 



12| . In the general quantum case, the additional equations determine the time evolution 
of the energy distribution function which in general can not be expressed in terms of the 
equal-time Wigner function. In this sense the BGR equations do not provide a complete set 
of equal-time kinetic equations. 

As we will discuss in this paper, this incompleteness has a more general aspect. As just 



mentioned, the equal-time Wigner operator is related to the covariant one by [12 



W(x,p) = I dEW(x,p), (1.1) 

where we wrote p = (E,p), E independent of p. As such it is only the lowest member of an 
infinite hierarchy of energy moments of the covariant Wigner operator: 

W J (x,p)= JdEE j W(x,p), j = 0,l,2,---, (1.2) 

with W (x,p) = W(x,p). Therefore, to set up a complete equal-time transport theory 
which contains the same amount of information as the covariant theory one needs dynamical 
equations for all the energy moments. Any covariant kinetic equation will thus correspond 
to an infinite hierarchy of coupled kinetic equations for its energy moments, i.e. for the 
equal-time Wigner operators Wj(x,p). 

This infinite hierarchy only exists for genuine quantum problems where the energy can 
exhibit quantum fluctuations. In the classical limit, the covariant Wigner operator satisfies 
the mass-shell constraint p 2 = E 2 — p 2 = m 2 , and the energy dependence of the covariant 



Wigner function thus degenerates to two delta-functions at E = ±£L = ±^m 2 + p 2 . The 



equal-time Wigner operator W(x,p) then splits into a positive and a negative frequency 
component, 

W(x,p) = W + (x,p) + W'(x,p) , (ft->0) (1.3) 

and all energy moments can be expressed algebraically in terms of these as 

W J (x,p) = W;(x,p) + Wr( x ,p) 

= E{,W + (x,p) + (-E p yW-(x,p), j = 0,l,2,... (fc-O). (1.4) 

The solution of the equal-time kinetic equations for W(x,p) thus also determines the dy- 
namics of all higher energy moments. Thus, in the classical limit, a simple zeroth order 
energy average of the covariant kinetic theory yields a complete equal-time kinetic theory. 

In the general quantum case, the higher order energy moments Wj(x,p), j > 1, contain 
genuine additional information and can no longer be expressed algebraically through the 
equal-time Wigner operator W(x,p). This means that in principle in the equal-time formu- 
lation we are stuck with the problem of solving an infinite hierarchy of coupled equations. 
Actually, there are two such hierarchies, one resulting from the covariant transport equa- 
tion ("transport hierarchy"), the other arising from the generalized mass-shell constraint 
("constraint hierarchy"). In practice this raises the problem of truncating the hierarchy in 
a physically sensible way. Since only the low-order energy moments of the covariant Wigner 
function have an intuitive physical interpretation, it turns out that physics itself suggests an 
appropriate truncation scheme. We will show that the hierarchies of moment equations are 
structured in such a way that the first few low-order moments form a finite and closed sub- 
group of equations which can be solved as an initial value problem, and that (surprisingly) 
all the higher order moments can be derived from these low-order moments recursively using 
only the constraint hierarchy, i.e. without solving any additional equations of motion. The 
equations from the transport hierarchy for the higher order moments are redundant. 

We first discuss on a general basis, starting from the covariant approach, the deriva- 
tion and truncation of equal-time hierarchies of kinetic equations. For illustration we then 
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consider in full generality the case of a transport theory for scalar fields with arbitrary 
scalar potentials. For this case everything can be worked out explicitly to arbitrary order 
of the moments. We give the subgroup of equations which fully characterize the first few 
low-order moments, prove the independence and redundancy of the transport equations for 
the higher order moments outside this subgroup, and obtain from the constraint hierarchy 
explicit expressions for all the higher order moments in terms of the solutions of the low- 
order subgroup. We then apply the general formalism to scalar and spinor QED. Here the 
equations have a more complicated structure, and we restrict our attention to the closed 
subgroup of equations for the lowest order moments, discussing the redundancy of the trans- 
port equations for the higher order moments and their recursive determination through the 
constraint hierarchy only for the first moments outside the closed subgroup of equations 
for the low-order moments. We will compare our results with the kinetic equations for the 
equal-time distribution functions obtained previously in Refs. [f7|,p|,[12| . Our final result will 
be a complete set of kinetic equations which can be implemented numerically as an initial 
value problem. 



II. GENERAL FORMALISM 

The 4-dimensional Wigner transform of the equation of motion for the covariant density 
operator leads to a complex, Lorentz covariant kinetic equation for the Wigner operator. It 
couples the one-body Wigner operator to two-body correlations Q, which in turn satisfy 
an equation which couples them to three-body terms, and so on. After taking an ensemble 
average this generates the so-called BBGKY hierarchy |13[ for the n-body Wigner func- 



tions. A popular way to get a closed kinetic equation for the one-body Wigner function 
(i.e. the ensemble average of the one-body Wigner operator) is to truncate the BBGKY hi- 
erarchy at the one-body level, by factorizing the two-body Wigner functions in the Hartree 
approximation. So far most applications of quantum transport theory have employed this 
approximation, and in the following we will also restrict ourselves to it. For us the mean 
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field approximation provides a crucial simplification, and at present it is not obvious to us 
how to generalize our results in order to include correlations and collision terms. 

For a scalar field in mean field approximation the complex equation for the self-adjoint 



scalar Wigner function can be separated into two independent real equations fl2f : 

G(x,p)W(x,p) = 0; (2.1a) 
F(x,p)W{x,p) = 0. (2.1b) 

The first equation corresponds to a generalized Vlasov equation; after performing the en- 
ergy average it generates a hierarchy of transport equations for the energy moments Wj (x, p) 
("transport hierarchy"). The second equation is a generalized mass-shell constraint; it gen- 
erates a hierarchy of non-dynamic constraint equations ( "constraint hierarchy" ) . Equations 
with the structure given in ( |2.1|) will be the starting point for our discussion of scalar field 
theories in Sees. |TXT| and |IV A| . Factors of p M in the dynamical operators G(x,p) and F(x,p) 



arise from the Wigner transformation of the partial derivative in the Klein-Gordon equa- 
tion. Since the latter contains at most second order time derivatives, at most two powers 
of po occur. In fact, G(x,p) is linear in p while F(x,p) is quadratic in p . This will be 
important below (see Sec. [11 B| ). 

For spinor fields the covariant Wigner function is a 4 x 4 matrix in spinor space which 
for a physical interpretation must be decomposed into its 16 spinor components W s . In this 
way the complex kinetic equation for the Wigner function matrix is split into 32 independent 
real equations for the self-adjoint spinor components |6|JT2|. These equations can be further 
divided into two subgroups according to their anticipated structure after performing the 
energy average: 



Hi 



G ss ' (x, p) W s ' (x, p) = , (2.2a) 

s'=l 
16 

£ F ss '(x,p) W s '(x,p) = , (s = 1, 2, . . . , 16). (2.2b) 



s'=l 



The first subgroup leads to equations containing only first order time derivatives and thus 
generates 16 hierarchies of transport equations for the energy moments of the spinor compo- 
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nents ("transport hierarchies").. The other subgroup which involves both first and second 
order time derivatives leads to a set of 16 hierarchies of constraint equations ("constraint 
hierarchies") for the equal-time moments of the spinor components. For the lowest energy 
moments, the spinor components W s (x, p) of the equal-time Wigner function, the details 
of this procedure were worked out in Ref. [O, and we will use these results in Sec. [IV C 



Since the original Dirac equation is linear in the time derivative, the dynamical operators 
G ss ' (x,p) and F ss '(x,p) contain at most single powers of po- in fact, the operators G ss '(x,p) 
are independent of po- 



A. Hierarchy of energy moments 

In this subsection we will concentrate for simplicity on a single covariant kinetic equation 
of the generic form 

G(x,p)W(x,p) = 0, (2.3) 

where G(x,p) contains at most two powers of po and of the space-time derivative operator 
<9 M , but an arbitrary number of derivatives with respect to the momentum space coordinates 
(see Appendix^. We will return to the full set of equations (|2.1|) resp. (|2.2|) in the following 
Sections. 

We begin by decomposing the energy dependence of the Wigner function into a basis of 
orthogonal polynomials hj(E): 

oo 

W(x,p) = J2^(x,P)h J (E) . (2.4) 

3=0 

The expansion coefficients Wj(x,p) are defined in the equal-time phase space. Using the 
orthonormality relation 

J d f i(E)h i (E)h j (E) = 5 ij , (2.5) 

where dfi(E) is the appropriate integration measure associated with the chosen set of poly- 
nomials hj, the equal-time components Wj(x,p) can be related to energy moments of the 
covariant Wigner function constructed with the basis functions hj(E): 
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Wj (x,p) = J dfM(E)h j (E)W(x,p) . (2.6) 

If the system has finite total energy the covariant spinor components must vanish in the 
limit E — > ±00. We will assume that they vanish at infinite energy faster than any power 
of E such that for any combination of integers i,j,m,n > 0, we have 

With exponential accuracy we may therefore restrict the energy integration to a finite interval 
—A < E < A. Introducing the scaled energy u = E/A we can thus use as our set of basis 
functions the Legendre polynomials 



h n {u) = ^±ip„( w ) (2.8) 

with the trivial measure dfi(uj) = du on the interval [—1,1]. 

As discussed above, the dynamical operator G(x,p) in ( |2.3p in general contains powers 
of E up to second order and an infinite number of energy derivatives d/dE. In terms of the 
new dimensionless energy variable we may thus write 

M 00 / o \n 

G(X, P )= E J2G™n(x,p)u m — , (2.9) 
m=0n=0 V^/ 

with M < 2. Substituting this double expansion into equation ( |2.3| ), multiplying by hi(ou) 
from the left and integrating over energy uj we obtain 

Y,G mn (x,p) C duh t (u)u m ^W(x,p)=0. (2.10) 

m,n 0uJ 

Inserting the expansion ( |2.4j) of the covariant Wigner function W(x, p, u), this can be written 



as 



where 



J2H l3 (x,p)w J (x,p) = 0, i = 0,l,2,..., (2.11) 

3=0 



M oo 

Hv(x,p) = EE^^.P) (2-12) 

m=0 n=0 

8 



with 



C™ n = f 1 duoh^^dlh^uo) . (2.13) 



It is easy to see from ( |2.13| ) that 



C™ = for n>j and i > j + m - n . (2.14) 

For j > i — m + n the coefficients are in general nonzero. Therefore, the sum over j in 
Eq. (12111) extends over the range j > max[0,i — M\. For each value of i, Eq. ( |2.11| ) thus 



contains an infinite number of terms. In this form Eqs. ( p.ll| ) are thus not practically 
useful. However, one can use the surface condition ( [2.7|) to rewrite Eqs. (|2.11| ) in such a 
way that each equation contains only a finite number of terms. Returning to Eq. ( |2.10| ) and 
integrating by parts, we can replace the integrand by 

h i (u)u m d:W(x,p)= (2.15) 

n-l 

£(-R (dlfaiu^d^w&p)) + (-rd2(h t (cu)cu m )w(x, P ) . 

1=0 



The contribution to the integral in ( [2.10|) from each term in the sum is fully canceled by 



the surface condition ( |2.7p , and only the last term in (|2.15|) survives. Inserting it into ([2.1 0|) 
and using again the expansion ( p.4| ) we find instead of Eqs. (|2.11| - pT3"D the following set of 
equations: 

i+M 

9ij(x,P)wj(x,P) = 0, i = 0,1,2,..., (2.16) 

3=0 



with 



M i+m 

g^p) = E E c^G mn (x,p) (2.17) 

m=0 n=0 



and the coefficients 



= i > /(2t + l)(2j + l) £ i duP j (cj)(-&r(P i (u,)u, m ) . (2.18) 

The latter can be determined recursively from = 5^ (which results from the orthogonality 
relation ( |2.5j )) by using the recursion relations for the Legendre polynomials, see Appendix 
Since the nonvanishing coefficients c™ n are now restricted to the domain 

9 



n < i + m and j < i + m — n , (2-19) 

the sum over j in (|2.16|) runs now only over the finite range < j < i + M. The first 
inequality in ( |2.19| ) was already used in ( |2.17| ) to limit the sum over n. 

The Pi{uS) are polynomials in uj of order i, and thus the equal-time components Wj(x,p) 
occurring in Eqs. fl2.16p are linear combinations of the energy moments Wk(x,p) of order 
k < i (see Eq. ( p.. 2D ). For each value of i, Eq. (|2.16|) thus provides a relation among 
the first % + M + 1 energy moments of the covariant Wigner function W(x,p) (including 
the zeroth order moment W(x,p) = v2wo(x,p)). As i is allowed to run over all positive 
integers, Eqs. ( |2.16|) form an infinite hierarchy of relations among the energy moments of the 



covariant Wigner function. Each covariant equation of the type (|2.3|) generates its own such 
hierarchy. Only the full set of these infinite hierarchies of moment equations constitutes a 
complete equal-time kinetic description of the system under study. 

B. Truncating the hierarchy 

In order to discuss possible truncation schemes we must return to the complete set of 
covariant kinetic equations. Let us concentrate here on the scalar case, Eqs. Q2.1|) , and write 
down the two resulting hierarchies of moment equations as 

i+M 

9ij{x,p)w j {x,p) =0, (2.20a) 

3=0 
i+M+1 

Y Ux,p)w j (x,p) = 0, (i = 0,1,2,...). (2.20b) 

3=0 



In writing down the upper limits of the sums we already used that F(x,p) in ( p. lap contains 



one power of po more than G(x,p) in ( |2.1a ). For the scalar field case one has M = 1. For 



the spinor case one obtains from ( |2.2| ) a similar set of equations with M = 0. 

Let us now try to truncate these hierarchies for the moments Wj at some order j max . 
The equations from the "transport hierarchy" (|2.20a|) with hierarchy index % < I t involve all 



moments Wj with < j < I t + M, i.e. the first I t + M + 1 moments (including the lowest 
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moment with index 0). Similarly, the equations from the "constraint hierarchy" ( |2.20b|) with 



hierarchy index i < I c involve the first I c + M + 2 moments < j < I c + M + 1. For a 
closed set of equations both hierarchies must be truncated at the same order j max , i.e. we 
must have 

I t + M = I c + M+ l= j max . (2.21) 

Truncating in this way we are left with Jj+1 equations from the transport hierarchy and J c +1 
equations from the constraint hierarchy. In order solve them the number of equations must 
at least equal the number of moments. However, if there are more equations than moments, 
the system may be overdetermined, and therefore we would like to require equality of the 
number of equations and moments: 

I t + 1 + I c + 1 = j max + 1 . (2.22) 

The two conditions ( |2.2ip and ( 2.22 ) have a unique solution: 



I t = I c + l = M, Jmax = 2M, (2.23) 

which yields M + 1 transport and M constraint equations for the first 2M + 1 energy 
moments. Smaller values of j max don't yield enough equations, and larger values lead to 
an (at least superficially) overdetermined system of equations. For spinor fields (M = 0) 
the truncated set involves only one transport and no constraint equation; for each spinor 
component it gives a single kinetic equation (the BGR equation 0) for its lowest energy 
moment, the equal-time Wigner function W s (x,p). For scalar fields (M = 1) the truncated 
set contains two transport equations and one constraint for the three lowest order moments 
w , Wi, w 2 . 

If we go beyond this minimal closed subset of moments and equations, we get two equa- 
tions for every additional moment, one from the transport hierarchy and one from the 
constraint hierarchy. As we will show explicitly in the next Section, in the constraint hier- 
archy (|2.20b| ) the highest moment always comes with a constant coefficient. As we increase 
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the hierarchy index % in Eqs. (|2.20|) , at each step the newly occurring moment can thus be 
explicitly expressed in terms of the already known lower order moments using the corre- 
sponding constraint equation from ( |2.20b| ) . As we will discuss, these higher order constraint 
equations contain important physics. But in addition, at each step there is also a dynamical 
equation of motion for the new moment from the transport hierarchy ( |2.20a| ). How can the 
two equations be consistent? The answer is that this transport equation is not an indepen- 
dent new equation, but (with some algebraic effort) can be expressed as a combination of the 
lower order equations which have already been used. Our proof of this fact uses explicitly 
the structure of the dynamical operators G mn and F mn . It involves cumbersome algebra, 
and only for scalar fields with only scalar potential or mean field interactions we have been 
able to find a general proof. For scalar and spinor QED the proof is still incomplete, and 
we will only demonstrate the first step for the 2M + 2 nd moment. A completion of the proof 
presumably requires a so far missing deeper insight into the general dynamic structure of 
the moment equations and their relation to the underlying covariant theory. 



III. SCALAR FIELD THEORY 

In this Section we will give an explicit and complete discussion of the moment hierarchy 
for the simplest case of a scalar field theory in Hartree approximation. We exemplify the 
truncation of the hierarchy and the recursive computation of the higher order moments 
beyond minimal truncation. The discussion in the following Section for the practically more 
relevant case of QED will be technically more involved and, unfortunately, also less complete. 

A. Covariant kinetic equations 

Consider the Klein-Gordon equation with a scalar potential U(x): 

(dl + ml + U{x)) 4>{x) = . (3.1) 
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The covariant Wigner function is the four-dimensional Wigner transform of the covariant 
density matrix g(x,y) = (g{x,y)): 

W(x, p) = J d A y e ipy g(x, y) = J d A y e ip ' y (x + §) 0+ (x - §) ) . (3.2) 

To derive the kinetic equations for the scalar Wigner function, we calculate the second-order 
derivatives Qc^ + d^j and (§c^ — c^) of the covariant density operator, and then employ 
the Klein-Gordon equation ( O ) and its adjoint. After taking the ensemble average and 
performing the Wigner transform we obtain two complex kinetic equations: 

(R 2 - P 2 + ™o + U (x - %d p ) - ip-d x ) W{x, p) = , (3.3a) 
-P 2 + ™o + + ify) + W-dx) W{x,p) = . (3.3b) 

Since the scalar Wigner function is real, adding and subtracting these two complex equations 
yields two real equations of the type ( J2.1[ ) . After reinstating % the corresponding operators 
G and F are given by 

G(x,p) = hp-d x + lmM 2 (x,p) , (3.4a) 

h 2 

F(x,p) = -p 2 + — d 2 x + ReM 2 (x,p) , (3.4b) 
where the mass operator M 2 is defined as 

M 2 (x,p) =m 2 + T, e (x,p) + i£, (x,p) , (3.5a) 

t, e (x,p) = cos (~Y^J U ( x ) > (3.5b) 

S (x,p) = sin f-^-J • ( 3 - 5c ) 

Here the triangle operator A is defined as A = d x -d p where the coordinate derivative d x 
acts only on the scalar potential U(x) while d p acts only on the Wigner function. 



B. Semiclassical expansion 



In the general quantum situation the particles have no definite mass due to quantum 
fluctuations around their classical mass shell and collision effects in the medium. In the 
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situation here with only an external potential this is illustrated by the mass operator M 2 . 
Only in the classical limit % — > it reduces to the quasiparticle mass 

ReM 2 (i,p) = m 2 (x) = m 2 + U{x) , (3.6a) 
lmM 2 {x,p)=0. (3.6b) 

In this case the constraint equation reduces to the on-shell condition 

(p 2 - m 2 {x)) W (x,p) = (3.7) 

for the classical covariant Wigner function Wo- The classical transport equation arises from 
the general transport equation at first order in %. The first order contribution to the mass 
operator is 

ReM 2 (x,p) = 0, 

1m M 2 (x,p) = hm(x) (d x m(x))-d p , (3-8) 

and we obtain the covariant Vlasov equation 

(p-d x + m(x)(d x m(x))-dp)W (x,p) = , (3.9) 

with a Vlasov force term induced by x-dependent effective mass term. For scalar fields there 
is no first order quantum correction to the operator F in ( |3.4| ), and from the zeroth order 
term we obtain the mass-shell condition for the first order Wigner function: 

(p 2 -m 2 (x)) Wi{x,p) = . (3.10) 

This discussion holds universally for arbitrary potentials. If, for instance, U(x) is generated 
by the scalar field (j)(x) itself in Hartree approximation, 

U( X ) = - C + \(4>(x)4> + (x)} = -C + X J ^YV(x : p), (3.11) 

with a mass parameter C and a coupling strength A, this model provides a useful tool for a 
dynamical description of spontaneous symmetry breaking |§. 
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C. The 3-dimensional dynamical operators 



We now perform the energy average of the covariant transport and constraint equations 
( j2.ip and construct the hierarchy (|2.20| ) of moment equations. The first step is the double 
expansion of the type ( |2.9|) for the covariant dynamical operators G(x,p) and F(x,p): 



G mn (x,p) 



1, n 







hAdt for m 

hp-'Vx — b Q for m = n = 

~h (S )" ^ d t a o) for m = 0, n ^ even 

-^(f)"^) form = 0,nodd 

else 



[x,p) 



-A 2 



for m = 2, n = 



^-<9 2 + p 2 + ml + cr e for m 



n 







(^e) 

-Mi) n (^o) 



for m = 0, n ^ even 
for m = 0,n odd 
else. 



Here 



(3.12a) 



(3.12b) 



a e (x,p) = cos (^-V^-VpJ U(x) , 
<7 (x,p) = sin (-V x -V p j U(x) 



(3.13a) 
(3.13b) 



are the three-dimensional analogies of the covariant operators S e and S Q in Eq. ( |3.5| ). Again, 
the spatial gradients act only on U(x), while the momentum gradients act on the equal-time 
Wigner functions (i.e. on the energy moments Wj(x,p)). 

The three-dimensional dynamical operators G mn (x,p) and F mn (x, p) must now be com- 
bined with the coefficients c™ n to obtain the dynamical operators §ij(x,p) and fij(x,p) 
which are needed in the transport and constraint hierarchies. This is done in Appendix [B]. 
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D. Minimal truncation 



The resulting transport hierarchy is truncated at It = M = 1, the constraint hierarchy 
at I c — M — 1 = 0. This yields the following equations for w (x,p), Wi(x,p), and w 2 (x, p): 



goowo + goiwi = , (3.14a) 
9ww + gnw 1 + g 12 w 2 = , (3.14b) 
/oo^o + /oi^i + fo2W 2 = . (3.14c) 

The dynamical operators and are given in Appendix and Eqs. (|3.12p . Reexpressing 
Wj in terms of the energy moments Wj from Eq. ( |1.2| ), 



w (x,p) = j^=W(x,p) , (3.15a) 
Wl (x,p) = ^^Wi(x,p), (3.15b) 
w 2 (x,p) = (^^ 2 (x,p) - W(x,p)} , (3.15c) 



Eqs. Q3.14 ) can be rewritten as 



dtW^x.p) = - (p-V x - -a {x,p)j W(x,p) , (3.16a) 
d t W 2 {x,p) = - (p-V x - ^& {x,p)j W!{x,p) + ~ (d t a e (x,p)) W(x,p) , (3.16b) 
W 2 (x,p)= (d 2 t -Vl) + p 2 + m 2 + a e {x,p)^W{x,p) . (3.16c) 

Please note that all powers of the cutoff A cancel in the final expressions as they should. 

The two transport equations ( |3.16| a,b) do not decouple, not even in the classical limit 
h — > 0. To achieve decoupling one must return to the covariant equations in Sec. |111 A| and 



study their semiclassical limit as given in Sec. [llllj| before performing the energy average. 
Then the mass-shell condition ( |3.7| ) can be used to rewrite all higher order energy moments 
in terms of the zeroth order moment as explained in the Introduction, Eq. ( |1.4j ). With this 
information the constraint (|3.16cj ), in the limit h — > 0, becomes trivial, 
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W 2 (x,p) = El(x)W(x,p), E 2 p (x) =p 2 + m 2 (x), (3.17) 
while the two transport equations ( [3.16a| ) and (|3.16b| ) become identical and can be written 



in the form of a Vlasov equation for the charge density (see Sec. [TV]): 



d t Wi(x,p) + ( J--V, - V x E p ■ V p j W x {x,p) = . (3.18) 

The reason why the information contained in Eq. ( |1.4j ) cannot be easily recovered directly 
from the 3-dimensional transport and constraint equations is that in their derivation, through 
Eq. (p.!5| ), we made heavy use of partial integration with respect to the energy. In the 
classical limit this has the unfortunate effect of spreading the information contained in the 
on-shell condition over the whole infinite hierarchy of 3-dimensional constraint equations. 

Although we have always talked about Eq. ( |3. 16c|) as a "constraint equation" it is clear 
that, as far as solving the minimal subset ( |3.16|) of equal-time kinetic equations is concerned, 
this terminology is only adequate in the classical limit h — > 0. In general quantum situations 
it is a second order partial differential equation for the lowest order moment W(x, p) which 
must be solved as an initial value problem together with the first order partial differential 
equations (|3.16a|) and (|3.16b|). 



E. Higher order moment equations 

The next higher moment u>3 is determined by the third equation in the transport hier- 
archy and the second equation in the constraint hierarchy, 

#20^0 + #21™1 + #22^2 + <?23^3 = , (3.19a) 

fww + fuwt + fi 2 w 2 + /13W3 = . (3.19b) 
With the dynamical operators from Appendix || and Eqs. (|3.12|) and 

w 3 (x,p) = -L^| (jgWfap) - |Wi(rr,p)) (3.20) 
we obtain (using Eq. ( |3.16a| )) 
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d t W 3 (x,p) 



W 3 (x,p) 



P-^x - ^a D (x,p) ) W 2 (x,p) + (d t a e (x,p)) Wi(x,p) 



n 



-(d 2 a (x,p)) W(x,p) 



4 

'h 2 

4 ^ 



— ( d 2 - V 2 ) + p 2 + m 2 + a e (x, p) W^x, p) 



■~ (a t 5" (x,p)) W(a:,p) . 



(3.21a) 



(3.21b) 



By substituting ( |3.19b| ) into (|3.19a|) and taking into account the following commutators: 



[Gio, -Foo] — — 2Goi-^2o > 



[Goo, F o\ 

h 2 

[—d 2 , cr e/o ] = —da e/o -d + ^—d 2 a. 



02-^20 j 

h 2 



r e/o ) 



[p-V x , a e/o ] = p-V x a e/o + -V x a o/e -V x 

h 2 

[P 2 , v e /o\ = -hp-V x a o/e + —V 2 a e/o , 



(3.22a) 
(3.22b) 
(3.22c) 
(3.22d) 

(3.22e) 



the third transport equation ( |3.19a| ) can be rewritten in terms of the first transport equation 
( glial) as 



foo (,900W0 +§01^1 ) = 0. 



(3.23) 



This implies that the transport equation ( |3.19a|) for w 3 is redundant. The third-order mo- 



ment W 3 is completely determined in terms of the solutions of the minimal subgroup ( 3.16 ) 
by the constraint equation (|3.21b| ). It arises from Eq. ( |3.19b|) by noting that /13 is a constant, 



/13 = C 3 — — |^|A 2 (see Appendix |B|) , and solving for w 3 : 



w 3 



~T (flQWo + fllWi + f 12 w 2 ) ■ 



(3.24) 



Note that, in contrast to (|3.16cj ), Eq. (|3.21b|) does not require solving a partial differential 
equation because everything on the r.h.s. is known from the solution of Eqs. ( |3.16| ). 

The above procedure can be extended to all higher order moment equations, by repeat- 
edly using the commutators listed in ( |3.22| ). In general one finds that a transport equation 



E 9ij W j = 
3=0 



(3.25) 



with i > 2 can be re-expressed in terms of the first (i — 1) transport equations as 

i-2 fj+i \ 

EU- E^* =°- ( 3 - 26 ) 

j=0 \fc=0 / 

Thus, except for the first two, all transport equations are redundant. The higher order 
moments W{ with i > 2 can be computed from their constraint equations 

1 . 

™< = (3.27) 

with the constant C{ given by 

C« = fi-2,i — f ' = • (3.28) 

(2« — 1)^(2* — 3)(2i + 1) 



IV. APPLICATION TO QED 
A. Scalar QED 

In this Section we discuss the application of the general formalism developed in Sec. |I|to 
QED. Since some of the equations will be rather lengthy we will economize on the notation 
by dropping all factors of H. The latter are correctly given in Refs. JT2]JT7[] to which we refer 
in case of need. 



We begin with the case of scalar QED with external electromagnetic fields. In Ref. [12 
we discussed the semiclassical transport equations for this theory by energy averaging the 
semiclassical limit of the covariant transport equations. In this subsection we will derive the 
general equal-time quantum transport equations by performing the energy average without 
any approximations. In the following subsection the result will be compared with the cor- 
responding equations derived by directly Wigner-transforming the equations of motion for 
the equal-time density matrix. 

In scalar QED the scalar field obeys the Klein-Gordon equation 

(D^D 11 + m 2 )(j)(x) = . (4.1) 
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From the corresponding Lagrangian density the conserved current and canonical energy- 
momentum tensor are derived as 

j^x) = ie {{D^{x)) <j>\x) - 4>(x)Dl4>\x)) , (4.2a) 
T^x) = {D^{x)){DW{x)) + {D u <f>(x))(D^{x)) - y^d-d (<P(x)<P\x)) , (4.2b) 

where = + ieA^(x) is the covariant derivative and its adjoint. 

Trying to express the components of the above energy-momentum tensor as moments of 
the covariant Wigner function results in rather complicated expressions. Simpler ones will 
be obtained after first subtracting the following total derivative terms 



\d»d v (<j>{x)<j>\x)) - y^d-d ((f>(x)tf(x)) 



(4.3) 



and employing the identity 

d$ v (0(x)0 f (x)) = (4.4) 
(D^(x))(DW(x)) + (D u <j>(x))(Dl^(x)) + (D,D u( j>(x))^(x) + ^(Ml^x)) 



to get a new tensor 



(4.5) 



| ({D^{x)){Dl<j>\x)) + (D„<f>{x)){DW(x)) - (D^ixWix) - 0(x)(Dt jD t t (x)) ) . 



Following an analogous procedure as in Sec. Ill A| (see [12] for details) one derives two 
covariant kinetic equations of type (|2.1|) for the covariant Wigner function 

,1/2 



W(x,p) = J d 4 ye ip ' y ^0 (x + §) exp ie J ^ ds A(x + sy)-y 4> ] (x 
The corresponding covariant dynamical operators G and F are given by 
G(x,p) = fl^(x,p)b^x,p) , 

F(x,p) = ^-D ll {x : p)b^{x : p) - W{x,p)%(x,p) + m 2 , 
fl M (x, p)=p fJi -ie / ^ ds s F^ v (x - isd p ) d v p , 

"2 

i 

D^{x, p)=d ll -e j^ds F^(x - isd p ) d u p . 

"2 
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(4.6) 

(4.7a) 
(4.7b) 

(4.7c) 
(4.7d) 



pvv _ qua" _ Q v A? i s the electromagnetic field tensor. 

The structure of the equal-time transport theory for scalar QED is very similar to that 
for a scalar potential U(x) which we considered in the previous section. The difference 
resides solely in the expressions for the dynamical operators G mn and F mn . In Appendix |C| 
we provide the double expansions ( |2.9| ) for the basic operators IP and -D M as well as for G 
and F. 

Upon expressing the field product (jxft and its covariant derivatives in terms of the 
covariant Wigner function (|4.6| ). 

4>{x)4>\x) = J -j?*W{x, P ) , (4.8a) 
(D^(xW(x) = J - W(x,p) , (4.8b) 

<j>(x)(D^(x)) = J -0- (\d, + ip») W(x,p) , (4.8c) 



the phase-space densities of the charge current j M and the energy-momentum tensor T^ v are 
simply given by 

f{x,p) = 2ep ti W(x,p) , (4.9a) 
T^(x,p) = 2p ti p v W(x,p) . (4.9b) 

The factors 2 in these expressions account for the contributions of particles and antiparticles. 
After performing the energy average these equations translate into relations between the first 
three moments w ,Wi, w 2 and the equal-time phase-space distributions for the scalar density 
W(x,p), the charge density p(x,p), and the energy density e(x, p): 

w (x,p) = J^fi W ( x >P) > (4.10a) 
wi(x,p) = \h^J5P( x >P) > (4.10b) 
w 2 (x,p) = 2^J~ (-^e(x,p) - W{x,pf) . (4.10c) 
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The charge current density j(x,p) and the momentum density P(x,p) can be expressed in 
terms of W and p as 

P 

j(x,p) = 2epW(x,p), P(x,p) = -p(x,p) . (4.11) 



The subgroup fl3.14p determining the moments wq, W\ and W2 can thus be equivalently 



rewritten as transport and constraint equations for W, p and e: 

-dtp + 2(rf> + Tr-d)W = , (4.12a) 

d t e + -(d> + iz d + D)p + 2{fr D + tt-I - d t A + G-d)W = , (4.12b) 

e 

' J2 U 2 \ (A2 ^-2\ . _2 1 J A tit 27T 



e = 2 f -(rf t 2 - d ) - (tt 2 - 7T 2 ) + m 2 + A) W - — -p , (4.12c) 

The expressions of and were obtained via the relations given in Appendix || from the 
equal-time operators ( |C2| ) in Appendix y. We have also used the commutators 

[4 vr ] = D , (4.13a) 
[4 A] =2E . (4.13b) 

Due to the line integrals over s in the operators 11^ and which guarantee the gauge invari- 
ance of the formalism the equal-time operators G mn and F mn for QED are much more 
complicated than Eqs. (|3.12|) for the case of a scalar potential. This is the origin of the more 



complicated structure of the minimal subgroup ( |4.12| ) of 3-dimensional kinetic equations. 
Please note the sequence of the operators in Eqs. (f4.12|) : in particular the generalized time 



derivatives dt act on everything following them. Eqs. ( |4.12| ) are thus much more intricately 
coupled than the corresponding equations ( |3.16|) for the pure scalar case. 

It is instructive to integrate Eqs. ( |4.12| ) over p to obtain equations of motion for the 
corresponding space-time densities: 

d t p{x) + V-i(x) = , (4.14a) 
d t e(x) + V-P(x) - E(x)-j(x) = , (4.14b) 
e(x) = 2 J ^5 (j-(d? - V 2 ) +p 2 + W(x,p) . (4.14c) 
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The first two equations express the conservation of electric charge and of energy-momentum 
while the last equation gives the energy density in terms of the scalar equal-time phase-space 
density W(x, p) including quantum corrections. (The factor 2 again accounts for particles 
and antiparticles.) In Eq. ( |4.14b| ) P(x) = f(d 3 p/(2ir) 3 )P(x,p) is the momentum density in 
coordinate space, and the last term describes the conversion of field energy into mechanical 
energy by the work done by the electric field on moving charges. 

Let us now consider the next two equations in the hierarchy, the transport and constraint 
equations ( 3.19|) for the third-order moment W3. Again the constraint equation can be di- 



rectly solved for W3 in terms of the solutions Wq, Wi, W2 of Eqs. (|4.10|) , (|4.12|) (see Eq. (|3.24|) ) 



Our task is to show that the transport equation for W3 is redundant. To this end we sub- 
stitute the constraint ( |3.19b| ) into the transport equation Q3.19a| ) and use the scalar QED 
analogue of the commutators ( |3.22| ), namely ( |4.13|) and 



[4 B) = 3F, (4.15a) 

[4 tt] = -I + V x 7f , (4.15b) 

[J, tt] =-V x -G, (4.15c) 

[d, 7T ] = V X 7T , (4.15d) 

[4 d] = -eV x -V p / 2 i ds E(x + isV p , t) , (4.15e) 

"2 

[tt, 7r ] = -ie j^dss E(x + isV p , t) + eV x / * ds s 2 E(x + isV p , t)-V p , (4.15f) 

as well as the identity 

J\ds (2is + ^V X -V P - s 2 V x -V p ^j E(x + isV p ,t) =0 , (4.16) 

to rewrite ( |3.19a| ) in the form 

foo (Qoo^o + Qoiwi) + /01 (^10^0 + gnwx + g 12 w 2 ) = . (4.17) 

Note that in the derivation of the commutators ([4.15|) we used Maxwell's equation 

d^ u = , (4.18) 
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for the dual field strength tensor = ^e^ uljp F ap . The identity Q4.16Q (which has no 
analogue in a theory without gauge invariance) is proven by expanding the electric field 
E(x + is'Vpjt) = e ls ' Vx '' Vp E(x) and integrating term by term. Eq. ( |4.17|) expresses the 



transport equation for in terms of the transport equations for W\ and w<i- This proves 
that it is redundant. 

From the comparison of (|4.17| ) with (|3.23|) we observe that in the case of scalar QED the 



third transport equation is expressed in terms of the first and second transport equations, 
while only the first one occurs in the case of scalar potentials. This difference results from the 
line integrals which occur in the gauge theory. In a transport theory with gauge invariance, 
the kinetic momentum is not but fl M . The energy average of the zeroth component of the 
second term in ( |4.7c|) yields the equal-time operator tc given in (|C2|) , and n = — Ir/oi i n 



turn generates the coefficient in front of the second bracket in Q4.17 ). In a transport theory 



without gauge freedom all coefficients /j^+i vanish due to the absence of linear terms in p M 
in the covariant constraint equation. 

We have not had the patience to carry the above considerations to higher orders in the 
energy moments. The corresponding calculations become extremely messy. Based on the 
experience with pure scalar theory we expect all higher order transport equations to be 
redundant, but we have failed to discover a simple calculational technique which permits us 
to prove this in an elegant way. 

B. Comparison with the Feshbach-Villars approach 

In Ref. H a set of equal-time transport equations for scalar QED was derived by directly 
Wigner transforming the equations of motion for the (non-covariant) equal-time density 
matrix. Since the Klein-Gordon equation ( |4.1| ) contains a second-order time derivative, 
its direct translation into 3-dimensional phase space does not lead to a sensible transport 
equation which should have only first-order time derivatives. Therefore the 3-dimensional 
approach exploits the Feshbach-Villars representation []l4j of the field equations of motion 
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which contains only first-order time derivatives. The price to pay (in addition to the loss 
of manifest Lorentz covariance) is the introduction of an auxiliary field which results in a 
rather complicated 2x2 matrix structure of the scalar Wigner function. For the energy 
averaging method advocated here there is no such problem. In this case second-order time 
derivatives appear only in the constraint hierarchy, and the transport hierarchy contains 
only first-order time derivatives. 

In terms of a two-component Feshbach-Villars field 

$=^j, (4.19a) 

V, = i(0+-9 t 0-^V) , (4.19b) 
2 V m m J 

X = \U--M+ 6 -^A , (4.19c) 
2 V m m J 

we define the density matrix 

QF{x,v) = (*(x + ±t) e--C^ A ^y'%^ x -p^ . (4.20) 

Its three dimensional Fourier transform with respect to y yields the equal-time Wigner 
function in Feshbach-Villars representation: 

W F (x, p)=J d 3 y e~ ipy g F (x, y) . (4.21) 

The equations of motion for the Wigner function Wf are a direct consequence of the field 
equation 

idt$ = (J- (-iV - eA) 2 (a 3 + ia 2 ) + ma 3 + eA ) $ (4.22) 
\2m / 

and its adjoint. Here 02 and 03 are the well-known Pauli matrices. We calculate the second- 
order derivatives of the density matrix qf with respect to x and y: 

Qv x - Q F (03 - icr 2 ) 

= 2m U(x + f ,t) e -*/-S a *^"^.*>-» ^t (a . _ u t ) ^ % +ma 3 + eA (x - f ,t) 



25 



-lie J ^^ds (~ - s) yxB(x + sj/,t) (-V^ - Vj,) £f (cr 3 - io 2 ) 
+e 2 i^J l/2 ds (I ~ s ) y xB (x + sy,t)^j £f (o"3 - icr 2 ) 

+ie J ^ ^ds (± - s) y ■ (V x+sy xB(x + sy,t)) g F (a 3 -icr 2 ) , (4.23a) 



= -2mU(id t - ma 3 - eA {x + f , tj) ®{x + §, t)) e~ ie ^ ds A{x+sy ^ y ^( x - U t ) \ 

-2ie J ' / ds (§ + s) yxB(x + sy,t) (|V X + V y ) (cr 3 + icr 2 )^ 

+e 2 (j^ds (§ + s) yxB(x + sj/,t)J (<x 3 + icx 2 ) £f 
/•1/2 . .2 

+ie J ^ 2 (| + s) y ■ (V x+sy xB(x + sy, t)) (a 3 + ia 2 ) g F . (4.23b) 

Note that we used the Feshbach-Villars equations for $ and <&t as well as Maxwell's equations 
for E = — VA — ^ and B = V x A. Subtracting the two equations Q4.23 ) we get a closed 



equation of motion for the density matrix g F (x,y); after a Fourier transformation with 
respect to y we obtain the following transport equation for the Wigner function in terms of 
the operators dt, d and 7r used in the last subsection: 

2md t W F + i (^\dt - 7r 2 ^ (W f (cr 3 - ia 2 ) - (cr 3 + icr 2 ) Wf) 

+7T • d(]V F (a 3 - ia 2 ) + (<x 3 + ia 2 )W F ) - 2im 2 {W F a 3 - a 3 W F ) = . (4.24) 

Expressing the equal-time Wigner function W F in terms of the Feshbach-Villars spinors || , 

W F = §(/ 3 + /atT! + /ia 2 + / <7 3 ) , (4.25) 

we recover the transport equations for the spinor components fi(x,p)(i = 0,1,2,3) which 
were first derived by Best, Gornicki and Greiner ||: 

mdtfo = -ir-d(f 2 + / 3 ) , (4.26a) 

mdtfr = - [\d 2 - tt 2 ) (f 2 + / 3 ) + 2m 2 / 2 , (4.26b) 

md t f 2 = (\d 2 - 7T 2 ) h + 7r-d/ - 2m 2 /! , (4.26c) 

mdth = ~ {\d 2 - 7T 2 ) h - n-dfo . (4.26d) 
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The current j^{x) from ( 4.2a|) and energy- momentum tensor T lxv {x) from ( |4.5| ) can be ex- 



pressed in terms of the Feshbach-Villars fields $ and <&t by using the transformation (|4. 19|) 



Similar to the equations (|4 . 8|) , we have relations between the field products their 



derivatives, and the Wigner function W F : 

$(x)$t(.x) = J -0-W F (x,p) , (4.27a) 

-iV-eA)$(xj)&(x) = J -0^(- l -V+p^W F (x,p) , (4.27b) 

*(x)((iV-eA)&(x)) = J^^V + p^Wpfap) . ( 4 - 27c ) 



Inserting these together with the spinor decomposition (|4.25|) into the expressions for and 



Tjj;, and comparing the integrands of the momentum integrations, one is led to the following 
identifications for the equal-time phase-space densities p, e, j and P: 

p F (x, p) = 2emf (x, p) , (4.28a) 

e F (x,p) = [p 2 - (f 2 (x,p) + h(x,p)) + 2m 2 f 3 (x,p) , (4.28b) 

j F (x,p) = 2ep(f 2 (x,p) + f 3 (x,p)), (4.28c) 

P F {x, p) = 2mpf {x, p) . (4.28d) 



Note that here, unlike the spinor decomposition for spinor QED ||[7] where each component 
of the Wigner function has a definite physical meaning, there is no obvious physical inter- 
pretation for the component f\. We also point out that the above expressions, especially 
the one for the energy density, differ from those given in Ref . |B],[T5|1 . The index F was added 
in order to point out that these phase-space densities are not equal to those in Eqs. Q4.10D , 



( [4. 11| ) - they generally differ by a total momentum-space derivative, and only their respec- 
tive momentum space integrals (i.e. the corresponding space-time densities) are guaranteed 
to be identical. The explicit relation between the two sets of phase-space densities is given 
in Appendix |TJ, Eqs. (pg|). 

Inserting the relations ( 4.28Q into the Feshbach-Villars equations of motions (|4.26[ ) one 



27 



finds, by choosing suitable linear combinations, the following equations of motion for the 
physical phase-space distributions pp and ep\ 

^d t pF + 2-k-dW = , (4.29a) 
1 . \ / \d 2 „ N 4 * , o V 2 



rf ; eF + .^..^ _ u_ _ + rf ;( p 2 _ w = ) (429b) 

e F = 2 (± $ - + 1 (tt 2 + p 2 - ^) + m 2 )^ . (4.29c) 

Here we used the shorthand W(x, p) for the combination 

W(x,p) = f 2 (x,p)+f 3 (x,p). (4.30) 

In Appendix it is shown that W(x,p) defined in this way indeed agrees with the scalar 
equal-time Wigner density (|4.10a|) of the Klein-Gordon field. In addition to the three equa- 
tions above, ( 4.26|) yields a fourth equation of motion for the unphysical Feshbach-Villars 



spinor component fx in terms of the scalar Wigner density W: 



f 1 = --Ld t W. (4.31) 
2m 



Apparently, the set of equations Q4.29 ) from the equal-time Feshbach-Villars approach 



and the set ( |4.12| ) derived in the energy-averaged covariant Klein-Gordon approach have 



a different structure. But in Appendix we show that both sets of equations are in fact 
equivalent. The structural difference is only due to the fact that the phase-space densities 
p,e in the covariant Klein-Gordon approach differ from pf,£f introduced via ( |4.28| ) in the 
Feshbach-Villars approach. This difference disappears in the classical limit where both sets 
of equations reduce to 

-dtp + 2p-dW = , (4.32a) 

dte + -p dp - 2ep-EW = , (4.32b) 

e 

e = 2(p 2 + m 2 )W. (4.32c) 

(Here we have used the classical operators d t = d t + eE-'Vp and d = V + eB x V p , without 
prefactors H.) For constant external electric fields they coincide even on the quantum level, 
yielding two decoupled ordinary differential equations in time for W and p: 
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d t p = , 

(d? + 4(p 2 + m 2 )^ + 4e£-p) W = . (4.33a) 

These equations were studied before in []5,|TT,l2"]). The first equation expresses charge con- 
servation in a homogeneous electric field while the second one is the well-known equation 
of motion for the charge current j OJT2I] which, according to the relation j = e , p --n, 
governs the time evolution of the particle density n by pair production in the electric field. 
This latter quantity thus comes out the same in both approaches. 



Finally, after integrating the Feshbach-Villars equations (|4.29|) over momentum one ob 



tains the same conservation laws ( [4.14D for the space-time densities as in the Klein-Gordon 
approach. The difference between the densities p,e,P and pf,£f,Pf in phase-space are 
not visible on the coordinate space level. Since all previous studies were done in one of 
the above limiting cases, the subtleties related to the exact definition of the phase-space 
densities p and e were apparently not noticed before. By accounting for it correctly both 
the equal-time Feshbach-Villars approach and the energy-averaged covariant Klein-Gordon 
approach are seen to be fully equivalent. 

Before ending this comparison between the two approaches we would, however, like to 
point out that dropping the derivative term in the definition (|4.28b ) of the energy density 



as done in Ref. [§] leads to wrong conservation laws after momentum integration; instead of 
( 14.141 ) one then finds 

d t p{x) + V-j'(z) = 0, (4.34a) 
d t e(x) + V-P(x) - E(x)-j(x) = —^-dt J tAeW(x :P ) , (4.34b) 



e{x) = 2 



/ (f ^ - +p 2 + m2 ) w ^ > ( 4 - 34c ) 



which contains an unphysical term on the r.h.s. of the energy conservation law and a spurious 
factor | in front of the Laplace operator in the definition of the energy density which breaks 
Lorentz covariance. 
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C. Spinor QED 



The case of spinor QED has been discussed in the context of the energy averaging method 
in Refs. |12|J17| , [19| . The discussion presented in those papers has, however, focused entirely 
on the equal-time kinetic equations for the lowest energy moment of the covariant Wigner 
function. Here we will reformulate the problem in terms of the equal-time hierarchy of 
energy moments as introduced in Sec. and also discuss the equations of motion for the 
higher order moments. 

Let us briefly review the relevant technical steps. We begin by performing a spinor 
decomposition || of the covariant Wigner function, separating in a second step explicitly 
the temporal and spatial parts J7| of the covariant spinor components: 

(4.35a) 



W(x,p) = - 



F{x,p) + ij 5 P(x,p) + ^V(x,p) + JsysAPfap) + ^a^S^(x,p) 



t [70^0(2;, p) +^ F 1 (x,p) + ij 5 F 2 (x,p) + F 3 (x,p) 

-l5l-G (x,p) - 7-Gi(x,p) + iry 5 j-G 2 (x,p) + >y 5 >yoJ-G 3 (x,p)} . (4.35b) 



Inserting the decomposition (|4.35a|) into the covariant (VGE) equation of motion || for the 
Wigner function W(x,p), 

(V(rt>,p) + \b^{x,p)) - m) W(x,p) = , (4.36) 

with fl^(x,p) and D^x^p) from Eqs. (|4. 7| c,d) , and separating real and imaginary parts one 
arrives at two groups of coupled covariant kinetic equations of type Q2.2|) , see Eqs. (74,75) 
of Ref. |jl2| . Performing the energy average then leads to two groups of equal-time kinetic 
hierarchies (with 16 such hierarchies of equations in each group) for the energy moments of 
the 16 covariant spinor components. 

Since for spinor QED M — 0, minimal truncation of these 2 x 16 hierarchies according to 
Sec. [11 B| results in one transport equation from each of the 16 transport hierarchies and no 



constraint equations. The minimal subgroup of equal-time kinetic equations thus consists 
only of the 16 transport equations for the 16 zeroth-order energy moments fi{x,p) and 
gi(x,p) of the covariant spinor components Fi(x,p) and Gi(x,p) {i = 0, 1,2,3): 
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dtfo + d-9i = , 




(4.37a) 


dth + d-9o + 2mf 2 = 


o, 


(4.37b) 


d t f 2 + 2ir-g 3 - 2m/i = 


= , 


(4.37c) 


dtf 3 ~ 2n-g 2 = , 




(4.37d) 


d t 9o + df x - 2tt x g x 


= 0, 


(4.37e) 


dt9\ + df - 2fr x g 


+ 2mg 2 = , 


(4.37f) 


d t g 2 + dx g 3 + 2-kf 3 


- 2mg 1 = , 


(4.37g) 


d t g3 ~dxg 2 - 2irf 2 


= . 


(4.37h) 



These 16 equations are identical with the BGR equations derived by Bialynicki-Birula, 
Gornicki and Rafelski @ by Wigner transforming the equations of motion for the equal- 
time density matrix. They determine the dynamics of the zeroth-order energy moments. 
The three-dimensional dynamical operators occurring in these equations are identical with 
the ones arising in scalar QED and are given in Appendix 0. 

The first-order moments fl(x,p) and gj(x,p) satisfy 16 transport equations derived [12 



from the first energy moment of the covariant transport equations ( 2.2a ), 



±=d t fl + -^d-g\ + bf + Ig x = , (4.38a) 

±J t f{ + ^=d-gl + A m /i + bh + I g = , (4.38b) 

\dtfl + ^f-9\ ~ A=mfl + Df 2 + 2Gg 3 = , (4.38c) 

■^dJs ~ 4=*-0a + Df 3 - 2G g 2 = , (4.38d) 

±=d t gl + ±=dfl + ^=Ttxg{ + Dg +J/ 1 -2Gx fll = 0, (4.38e) 

±=d t g\ + -j=dft -l=T T xg 1 + ^jfig\ + D 9l + If - 2G x g = , (4.38f ) 

±=d t g\ + -Ldxgl+ ^=kfl - ^=mg\ + bg 2 + ixg 3 + 2Gf 3 = , (4.38g) 

±J*g\ -±=dxg\- A^/i + bg 3 - I x g 2 - 2Gf 2 = , (4.38h) 



and 16 constraint equations derived |12fl from the zeroth energy moment of the covariant 
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constraint equation ( |2.2b|) : 



-f 1 

J_ i 
V3 9 ° 

J_ 1 

V3 92 



TT-f/l - 7To/o + mf 3 , 


(4.39a) 


7T-£o - fr f 1 , 


(4.39b) 


~^d-g 3 - 71-0/2 , 


(4.39c) 


^d-g 2 - 7r / 3 + mfo , 


(4.39d) 


-dxg 1 + tt/i - 7r gr + mg z , 
2 


(4.39e) 


-d X g + 7r/ - 7T C/i , 


(4.39f) 


~\ d h + 7T x g 3 - ir g 2 , 


(4.39g) 


^d/ 2 -tt x g 2 - 7r gr 3 + mgf • 


(4.39h) 



A discussion similar to that in scalar QED reveals [19|] that the transport equations (|4.38|) 



are not independent of the BGR equations (^4.37|) and the constraint equations ( |4.39| ). For 



instance, using Eqs. (|4.39| ), the transport equation for /q 1 can be expressed as an operator 
combination of the transport equations for / and g%: 

7T (d t fo + d-gi) - tt- (d t gi + df -2<frxg + 2mg 2 ) = . (4.40) 

Therefore, the first-order moments are fully determined in terms of the solutions of the BGR 
equations (f4.37|) for the zeroth-order moments by solving the constraint equations (|4.39|) . 



Again, we have not been able to find a simple proof that the same is generally true for 
all higher order energy moments, and we stopped here. We do, however, believe that such 



a proof must exist |2(J, and that therefore all higher order energy moments can be directly 
computed from the solutions of the BGR transport equations by solving the constraint 
hierarchy. 

It was shown in Ref. |]12| that in the classical limit the simple algebraic relation ( |1.4| ) 
changes the structure of the constraints ( f4.39| ) for the first-order moments and turns them 



into additional constraints for the zeroth-order moments (i.e. the equal-time Wigner func- 
tions). These extra constraints reduce the number of independent zeroth-order moments 
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from 16 in the quantum case to 4 in the classical limit. They are thus extremely important. 
As a result the BGR equations reduce to two decoupled Vlasov-type transport equations for 
the charge and spin distribution functions. In the general quantum case there are no such 
extra constraints on the equal-time Wigner functions One must solve all 16 coupled 



transport equations ( f4.37| ), but these solutions then fully determine also all higher order 
moments. These higher order moments have important physical meaning: the first-order 
moment of F (x,p), for instance, describes the energy distribution in phase-space. With the 
help of the constraint (|4.39a| ) it is given by 

e(x,p) =Tr J ^-E l0 W(x,p) = J^ti{x,p) 

^(mf 3 {x,p) -7r fo(x,p) +7C-gx(x,p)) . (4.41) 



V. CONCLUSIONS 

We have presented a universal method for the construction of equal-time quantum trans- 
port theories from the covariant quantum field equations of motion. It is based on energy 
averaging the covariant kinetic equations for the covariant Wigner operator (which is the 
Wigner transform of the covariant, "two-time" density matrix) and its energy moments. 
This procedure yields a hierarchy of coupled transport and constraint equations for the en- 
ergy moments of the covariant Wigner function, the so-called equal-time Wigner functions. 
We showed how, in the mean-field approximation, this hierarchy can be truncated at a rather 
low level, requiring the solution of only a small number of equal-time transport equations, 
and how the higher order energy moments (higher order equal-time Wigner functions) can 
be constructed from these solutions via constraints. 

The major advantage of the equal-time formulation of (quantum) transport theory is 
that the resulting transport equations can be solved as initial value problems, with bound- 
ary values for the equal-time Wigner functions at t = — oo which can be calculated from the 
fields at t = — oo. This is not the case for the covariant transport equations and the covari- 
ant Wigner function. The present paper thus provides an essential step in the direction of 
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practical computations of the dynamics of relativistic quantum field systems out of thermal 
equilibrium in the language of transport theory, i.e. in a phase-space oriented approach. The 
method presented here improves upon previous approaches by being much more systematic: 
we did not just focus on the lowest energy moments (which contain only a small fraction of 
the information contained in the covariant Wigner function), but we discussed and showed 
how to solve the complete hierarchy of moment equations. We had already before demon- 
strated for spinor QED that the non-covariant three-dimensional approach (which starts 
directly from the equal-time density matrix) yields an incomplete set of equal-time trans- 
port equations. In this paper we also discussed the case of scalar QED and showed that again 
the correct physical interpretation of the Feshbach-Villars spinor components f , . . . , f 3 in 
the direct non-covariant equal-time approach is not possible without a comparison to the 
energy-averaged covariant approach which we presented in Appendix |D|. We conclude that 
the only safe way of deriving a correct and complete set of equal-time quantum transport 
equations is by starting from the covariant formulation and taking energy moments of the 
covariant kinetic equations. The method can be generalized in a straightforward way to 
other types of interactions [^,|T^,|T^], including non-Abelian gauge interactions 

The structure of the hierarchy of equal-time quantum kinetic equations depends on the 
structure of the covariant field equations from which one starts. For scalar or vector theories 
with second order time derivatives one has to solve a coupled set of three equations for the 
three lowest energy moments, two resulting from the equal-time transport hierarchy and 
one stemming from the constraint hierarchy. For spinor theories with only first order time 
derivatives in the field equations one ends up with only one equal-time transport equation 
for the lowest energy moment of each spinor component of the Wigner function. All higher 
order energy moments can be determined from the solutions of these equations by solving 
constraints. 

Important further simplifications occur in the classical limit h — > 0: then all higher order 
energy moments can be expressed algebraically in terms of the zeroth energy moment, and 
the number of equations is drastically reduced. For scalar theories one obtains just one 
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Vlasov-type equation for the on-shell charge distribution function. For spinor theories one 
obtains two decoupled Vlasov-type equations (one scalar and one vector equation) for the 
on-shell charge and spin density distributions in phase-space. Again the only systematic 
way of deriving the constraints leading to these simplifications is by energy averaging the 
(classical limit of the) covariant transport equations. 

All results in this paper were derived in the mean field approximation, i.e. in the colli- 
sionless limit. It is generally known that including collision terms in the covariant transport 
equations leads to the appearance of non-localities in time ( "memory effects" ) in the equal- 



time transport equations pl| . It is not inconceivable that these memory effects lead to 
serious complications for the truncation of the equal-time transport hierarchy. This is cer- 
tainly an interesting and difficult problem for future studies. 
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APPENDIX A: THE COEFFICIENTS c"T 



From the definition of the coefficients c™ n in ( |2.18| ) and the orthogonality condition 
we see that 

c<% = . (Al) 



Furthermore, it is easy to see from ( [2.18|) that 



c™j n = for n > i + m or j > i + m — n . (A2) 



Therefore, if the finite sum over m has the upper limit M, the sum over j in Eq. fl2.16| ) 
extends only over the range j < i + M. 
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From the recursion relation for the Legendre polynomials 

uPi{u) = ^TT PmM + 27TT Pi - 1 ^ ) (A3) 

one easily derives the recursion relation 

m+l,n _ i + I m ,n . 1 m,n /a,\ 

C i,j ~ T, C i+l,j ^ / , , C i-l,j 



/ (2z + l)(2z + 3) v /(2i+l)(2i-l) 
for the coefficients c™ n . This allows to raise the first upper index m, starting from ( |A1| ). In 
particular we have 

~ 10 - ' (A5a) 



(2i - l)(2i + 1) 



= ~ ~ ^ +1 L_. = > (A5b) 



(2i + lty(2i-l)(2i + 3) 
which we will need in Appendix |B[ Similarly the second upper index can be raised by using 

the relations 

d u {P j+1 {u) - Pj-i(u;)) = (2j + 1)P» , i>-(±l) = (±l) j (j > 0) , (A6) 
which, for j > 2, lead to 



m.n+X I "I - 1 m,n+l . //n • i i\/o- i \ m i n / a rr\ 

c» j = W ^— 3 + V + l )W " X ) c m-i • ( A7 ) 

This expression is useless for j = and j = 1; these cases can be treated by another recursion 
relation which can be obtained by using ( |A6| ) on the Legendre polynomial with the index i: 



m,n+l /2z -(- 1 m ,n+l // • i 1 \ /o • i~N m.n / m— l,n /2z -|- 1 m —l,n \ i a o\ 

ClJ = V 27^3 ° l 2J ~~ V ( )( ~ } Ci lj ~ m I ~ V 2^3 Ci " 2j j ' ( } 
For i < 2 this must be used together with 

c ' ■ = V-l Cjjj.- (A9) 



-ij v H-l,j 

which follows from P~i(x) = Pj_i(s). 

The above recursion relations can be initialized with the following nonvanishing coeffi- 
cients for i,j, m,n < 1: 

JX) 00 _ 11 _ -. 10 _ 10 _ 1 „01__ /o li _ _ 9 /A in\ 

c 00 ~~ c ll ~~ c 00 ~~ 1 > c 01 ~~ c 10 ~~ ^/g > c 10 — v ° ) L ll ~~ z • l^- ±u J 
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APPENDIX B: THE OPERATORS g i3 AND fa FOR SCALAR THEORIES 

The three-dimensional dynamical operators §ij(x,p) and fij(x,p) needed in Eqs. fl3.14| ) 
and ( 3.19 ) are obtained from the definition ( |2.17| ) with the coefficients c"j n from Appendix 



Please note that for the upper limit in ( p. 171 ) for the sum over m is M = 1 while for 
it is M + 1 = 2. One finds: 

9oo — Gqo ~ Gn j (Bla 

0oi = 7/|^io > ( Blb 

&o = v/3 f-Goi + ~<?io + 2G 12 ) , (Blc 

£n = Goo - 2Gn , (Bid 

2 - 

#12 = -7^10 , (Ble 

§20 

g 2 i = Vl5(-Goi + ^-Gio + 3Gi 2 ) , (Big 



'15 

VE(3d 02 -d 11 -9d 13 ) , (Bif 



15 

§22 = Goo - 3G U , (Blh 
3 

#23 = "T^^IO , (Bli 



x 35 

foo = Foo — Fu + 2^20 + 2-F 22 , (Blj 

foi = (Ao - 2F 2 i) , (Blk 

2 - 

/02 = ' ( B11 

Ao = v/3 (-Ai + 3 Ao + 2F 12 - Ai - 6F 23 ) , (Blm 

fn = Foo ~ 2A1 + -Ao + 6A2 , (Bin 
5 



/ 12 = F 10 -3F 21 , (Bio 

vlo v ' 



/ 13 = ^Ao • (Blp 

The dynamical operators G mn and F mn are given in equations ( |3.12| ) for scalar mean field 
interactions and in Appendix |C] for scalar QED. Please note that in both cases F 2 k = for 
k 7^ and thus some of the expressions for above simplify. 
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We also note that in 



2 i—l+m 

/i-l.i+l = X! C r™,i+l-^nn (B2) 

m=0 n=0 



the only term which doesn't violate the second inequality in ( |A2| ) is the one with m = 2, n 
0: 



/i-i,t+i - m F 20 - , + J F 20 . (B3) 

(2i+lty(2i- l)(2i + 3) 



In the last equality we used ( |A5b| ). From Eqs. ( |3.12b ) or ( |CJ4g ) we thus see that /i-i^+i is 



just a constant. Following the same reasoning one also obtains 

1 i+m ^ i ~ 

9i-l,i = }2 22 C T-l,iGmn = Cj_! jGio = , = G 10 , (B4) 

m=0n=0 V /(2Z-1)(2Z + 1) 



where according to Eqs. fl3.12a|) or ( |C3d| ) G w is proportional to the time derivative d t resp. 

APPENDIX C: THE OPERATORS G mn AND F mn FOR SCALAR QED 

The elementary operators IT M and Z) M are extensions of the covariant momentum and 
the covariant derivative d^, respectively. Their double expansions in u and d/du are 

fr a A d B d 2 C d 3 _ . 

n ° = A " + " o + A^-AW-Aw + --- ' (cla) 

5o -*"Aa;"w + w + '" ' (Clb) 

- G d H d 2 . 

n = 7r -^rirT + ^ir3 + ---> ( Clc ) 



7" r5 / <9 2 

^ = -rf+T7TT + ^7^ + --- , (Cld) 



A du A 2 du 2 
id J d 2 
A&J + A^dcu 2 

with the equal-time dynamical operators, 
, 1 

ttq(x, p) = ieh / ds s E(x + is^V p , t)-V p , (C2a) 
7r(x, p) = p — ieh j ds s B(x + ish'Vp, t) x V p , (C2b) 
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d t (x,p 
d(x,p 
A(x,p 
B(x,p 
C{x,p 
D(x,p 
E(x,p 
F(x,p 
G{x,p 
H(x,p 



I(x,p 
J{x,p 



hd t + eh / ^ ds E(x + ish'Vp,t)-'Vp, 



hV x + e% i ds B(x + ishVp, t) x V p , 



d 



eh 2 dss 2 —E(x + ishV p ,t)-'V p , 
• 1 at ' 



ieh 3 - 1 



d 2 

ds s 3 ^—E(x + ishV pi t)-V pl 
i at 1 



eh 4 



2 
1 

,2 '' 2 



— J x ds s A —E(x + ishV p ,t)-V p , 
d 

ieh 2 J" ds s—E(x + ish'Vp,t)-'Vp, 
J-- at 

eh 3 rh , o d 2 . . _ . _ 

— J_ i dss 2 —E{x + tshVp,t)-V 1 



ieh 



4 „1 



at 2 ' 
a 3 



6 7-i at 6 



,2 /' 2 „2 ^ 



2 

.3 i 



2 
1 

-2 ''2 



a 2 

(is s 3 — i?B(x + ishV p , t) x V„ 



<9 



efe | i dss 2 —E(x + ish'Vp } t), 







eh 3 - l 



ieh J i ds s—B(x + isk'Vp, t) x'Vp — eh J i dsE(x + isKV p , t) 
d 2 

ds s 2 —B(x + ishV p , t) x V p 

+ieh 2 ds s—E(x + ish'VjjA). 
J-l dt p 



2 7-i 



(C2c) 
(C2d) 
(C2e) 
(C2f) 
(C2g) 
(C2h) 
(C2i) 
(C2j) 



eh I i ds s 2 —B(x + ishV p , t) x V p + ieh j i ds sE(x + zsftV p , t), (C2k) 
ze/i 3 



(C21) 
(C2m) 



(C2n) 



By substituting these expansions into the expressions of G and F in (|4.7|) and using the 
definition of G mn ( |2.9|) , we list some low-order G mn and F mn which will be used in the 
derivation of the kinetic equations for low-order energy moments, 



G o = Kod t + 7r-d, 

d 01 = -— Ck D + TT-I - Ad t + G d 
G02 = -r* (H-d + G I - 7T- J - Bd t -A£>- txoE 
G\o = A.d t , 



(C3a) 
(C3b) 
(C3c) 
(C3d) 
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G 
G 

G 



11 — 


-D, 




E 


12 — 


~A' 




F 


13 = 


A 2 "' 



and 








-^00 


^V^t tx ) 


f j4 + TT 2 — 7T" 2 ") + TIT 2 
V / ' 5 




-?01 = 




-{d t ,D})-2n A-{*, 


6} + 2B) , 


Fq2 = 


h ^ - 


f + K j}-{4e»- 


- (i 2 - G 2 - 


Fio = 


-2A7TQ, 






F u = 


-2A, 






F\2 = 


4 






F20 = 


-A 2 . 







(C3e) 
(C3f) 

(C3g) 



(C4a) 
(C4b) 
(C4c) 
(C4d) 
(C4e) 
(C4f) 
(C4g) 



APPENDIX D: RELATION BETWEEN FV AND KG APPROACH 

In this Appendix we relate the Wigner function of the Feshbach-Villars (FV) approach 
to that of the Klein-Gordon formulation, both on the covariant and equal-time level. 

We rewrite the FV Wigner function Q4.21|) as an energy integral of the covariant FV 
Wigner density W F (x,p), 

W F (x,p) = J ^W F (x,p) =J^J d 4 ye^g F (x,y), (Dl) 

where the covariant FV density matrix p F {x,y) is defined as 

g F (x,y) + e ieI -^ dsA{x+sy>v & (x - |) . (D2) 



With the transformation (|4.19 ) the covariant FV 2x2 density matrix can be expressed in 
terms of the KG fields: 
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Q F (x,y) = e ie f-V* dsA{x+syyv x 

i fl + ±(D° + -D°) + ^D° + D° , 1 + — (z?9 + D°) \D° + D 

*\1-±(D°+D°) --^D°,D° , l-±(Dl-B°) + ^DlD 

^ m \ t — J m t — ' m \ + —J mr + 



(x + )0t( x „), (D3) 



where x± = x ± | and 



£>± = 9° ± zeA (:r±) . (D4) 



Using the following identities: 



d + e ieJ%dsA( X+ sy).y = fe + ^ ^ ^ + ^ ^ + ^ A je fj^s A {x+ sy>y ? 

and Fourier transforming with respect to y one obtains the covariant FV Wigner density in 
terms of the covariant KG Wigner function: 



Here Do(x,p) and flo(x,p) are the time components of the operators given in ( |4.7| c,d) . 
Performing the energy average we finally obtain a relation between the equal-time FV Wigner 
function and the energy moments of the KG Wigner function: 

W F (x,p) = hl + ax)W(x,p) - —a 2 d t (x,p)W(x,p) 
4 4m 

+^(^+V*,pW*,p)) 

[1 - cxi) ( \d\ (x, p)W(x, p) + -e(x, p) + -7Tq(x, p)p(x, p) 



4m 2 v v-"-/ 2 



e 



-A(x, p)W{x, p) + tt 2 (x, p)W(z, p) J . (D6) 

Here W(x, p), p(x, p) and e(x, p) are the scalar, charge, and energy density in phase-space as 
defined in Eqs. (|4.10|) for the Klein-Gordon field, and the operators 7to(x,p) and dt(x,p) are 
defined in Appendix y. Comparing Eq. (06) with the Feshbach-Villars spinor decomposition 



( 4.25| ) we find the following relations: 
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2emf = P + 2en W , (D7a) 

h = -±-d t W, (D7b) 

h + h = W, (D7c) 

m 2 (/ 2 - / 3 ) = (\d 2 t + tt 2 - i) W + |e + ivrop . (D7d) 

Inserting them into Eqs. (|4.28| ) we obtain the following relations between the FV phase-space 
densities (|4.28|) and the KG phase-space densities (14.10Q , (|4.11| ): 



p F = p + 2e7r W } (D8a) 
Of = j = 2epW , (D8b) 
e F = \e + -U p + [J (d 2 - V 2 ) + p 2 + m 2 + vr 2 - i] W , (D8c) 

P F = P + 2pTi W = -pp F . (D8d) 

e 

One easily checks that the momentum integrals of the FV and KG phase-space densities 
(i.e. the corresponding space-time densities) agree with each other. 
Inserting these relations into the equations of motion ( |4.29| ) we find 

-d t p + 2(d t 7t + 7r-d)W = , 

e 

d t e + ^ (c2t7ro + 7r-d) p + 2 ({(d 2 t - dt)d t + tt 2 ^ + rf> 2 - d t A + m 2 rf t + 27r-dvr ^ W = , 
e = 2 Q(d 2 - d) - (tt 2 - 7T 2 ) + m 2 + a) W - . (D9) 



The first and third equation are now seen to agree with the corresponding equations (|4.12|) 



from the energy-averaged covariant KG approach while the equation of motion for the energy 
density still looks different. However, by acting with tt on the first and with d t on the third 
equation and then combining them with the second equation, using the commutators ( 4.13 ) 



and Q4.15Q , one indeed also recovers Eq. ( |4.12b ). This shows that the Feshbach-Villars 



equations (|4.29| ) are equivalent to the minimal subgroup of equal-time kinetic equations 
( |4.12j ) in the energy-averaged covariant Klein-Gordon approach. 
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